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Abstract: 

We present the design of, and a first analysis of data from, the atmospheric seeing monitor at the Australia Telescope 
Compact Array (ATCA). The seeing monitor has been operational almost continuously since May 2004 and every 
lOmin delivers a measurement of the atmospheric phase stability at the observatory. Its measurements can be used 
by observers to help deciding whether it is worth carrying out observations at millimetre wavelengths or whether 
a longer-wavelength backup project should be observed. We present a statistical analysis of the data recorded 
since September 2004 to characterize the annual variations in atmospheric path length fluctuations. Our analysis 
shows that in terms of phase stability, nights in spring, summer, and autumn are as good as, or better than, days in 
winter. We also find that the data imply that the turbulence in the lower few hundred metres of the atmosphere is 
predominantly responsible for the atmospheric seeing. 

Keywords: methods: atmospheric effects, instrumentation: interferometers 



1 Introduction 

The dominant source of phase errors in interferometric ra- 
dio observations at frequencies above 5 GHz are fluctuations 
in the tropospheric water vapour content along the line of 
sight of the individual interferometer elements. The effect of 
the fluctuations scales linearly with frequency, making ob- 
servations at frequencies of tens of GHz (at wavelengths less 
than about one cm) particularly prone to atmospheric phase 
changes. For a radio interferometer operating at frequencies 
of up to 115 GHz such as the Australia Telescope Compact 
Array (ATCA) it is therefore desirable to identify the peri- 
ods within a year that are suitable for observations at these 
high frequencies. A well-tested means of characterizing the 
phase stability above a radio observatory are two-element in- 
terferometers observing unmodulated beacons of geostation- 
ary satellites (e.g. lRadford et alll99rllHiriart et all2002l) . 

1.1 Atmospheric seeing in the radio regime 



A exc ellent introduction to radio seeing can be found in Thompson 
1200 ll It is generally assumed that the water vapour distribu- 
tion in the troposphere can be described as follows. Water 
vapour does not mix well with air because the air temperature 
is close to the condensation point of water. Its distribution 
can be approximated as nested parcels of air with a variety 
of water vapour densities. Kinetic energy is transferred from 
larger parcels to smaller and smaller parcels until it is dissi- 
pated by viscous friction, a situation known as Kolmogorov 
turbulen ce. A sketch of the situation can be found in Fig- 
ure II in Masson 1 1994). The rms phase fluctuations (<7</>), 
measured over a long period, are expected to exhibit a power 
law relationship with baseline length, b: 



for the rms phase. When the layer is deep compared with the 
interferometer spacing, three-dimensional theory is appropri- 
ate, and a power law exponent of j3/2 = 5/6 is expected. 

Furthermore, it is assumed that the mixture of air and wa- 
ter vapour is "frozen" in the atmosphere and is transported 
relative to the ground at a particular speed, v s . This wind 
speed aloft couples the temporal and spatial behaviour of 
the phase screen, and v s therefore is called the phase screen 
speed. 



2 Design of the seeing monitor 

Here we report on the development of a seeing monitor at the 
ATCA, and present an analysis of some data. 

The ATCA seeing monitor is an interferometer on a 230 m 
east-west baseline and tracks the 30.48 GHz beacon on the 
geostationary communications satellite, OPTUS-B3, at an el- 
evation of 60° . Each element of the interferometer consists 
if-aJ.8 m prime-focus paraboloid equipped with an uncooled 
fc „ ■ nois 
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Generally the turbulence is considered to be confined to a 
specific layer. When the layer is much thinner than the inter- 
ferometer spacing, two-dimensional turbulence theory is ap- 
plicable, which suggests a power law exponent of/3/2 = 1/3 



9 M- noise amplifier (LNA). The FWHM of the reflectors' 
main lobes is 0.4°. This is large compared to the satellite's 
diurnal motion on the sky (0.1°), and hence the reflectors are 
stationary and do not need to track. 

The power arriving at the feed horns is approximately 
3 ■ 10~ 13 W, or -125 dBW. The signals are received and am- 
plified in the primary foci, then transported to a screened 
temperature-controlled enclosure next to the antenna, where 
the signals are down converted in two steps to a frequency of 
240 MHz. They are then sent to the ATCA's screened room 
via optical fibres, where they are mixed down to a frequency 
of 50 kHz, and analysed. 

The design of the system was challenged by the poor 
signal-to-noise ratio, and the drifting beacon frequency due 
to the satellite's radial diurnal motion: this motion makes the 
beacon frequency drift through approximately 5 kHz (peak- 
to-peak) each day. 

The problem of the changing frequency was addressed 
by coupling the local oscillator frequency to the satellite bea- 
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con. This is accomplished with a Stanford Research Sys- 
tems SR5 10 lock-in amplifier with a very narrow-band (2 Hz) 
notch filter, which tracks the satellite beacon frequency auto- 
matically. 

The signal-to-noise ratio was improved by generating a 
clean signal coupled to the notch-filtered frequency and using 
this signal as a reference to which the other antenna's signal 
is to be compared. This is realised with a Stanford Research 
Systems SR830 lock-in amplifier. The phase measurement is 
then carried out with a second SR5 10. 

In addition to the atmosphere, the interferometer phase 
is affected by two other phenomena: both the small motion 
of the satellite and thermal drifts in the receiving system in- 
troduce extra phase responses. FigureQshows a comparison 
between the measured phase and that predicted purely from 
the satellite's motion. The residuals after subtracting off this 
component clearly show both the atmospheric and thermal 
components. 

3 Observations and data analysis 

3.1 Operation of the seeing monitor 

The integration times of the seeing monitor's lock-in am- 
plifiers are set to 100 ms, however, the phase is read and 
archived by the observatory monitor system only every 5 s. 
The standard processing of the data are as follows: the differ- 
ence between two successive phase measurements is found, 
the standard deviation of this difference is computed and the 
measurement is converted into a path length using Eq. J3} 
before archiving. By taking the standard deviation (as dis- 
tinct from rms) of differences, any response to a linear drift 
within the averaging period is eliminated. Thus the effect 
of the satellite motion and thermal phase drifts are removed. 
This standard deviation is computed over 10 min of data. As 
lOmin is long compared to the wind-crossing time over the 
seeing monitor baseline (the crossing time is usually no more 
than one minute), this is a statistically adequate period for es- 
timating the rms phase. 

The real-time measurements and plots from the archive 
are accessible to staff and observers. As an example, the path 
length fluctuation from a seven day period are reproduced in 
Fig.H 

3.2 Our analysis 

The analysis presented in this paper is more detailed than the 
analysis carried out routinely at the observatory and yields 
more insights into the properties of the atmosphere. 

We have converted the raw seeing monitor phases into 
path length measurements in microns and divided them into 
sections of 30 min. A second-order polynomial was sub- 
tracted from the measurements to eliminate changes due to 
satellite motion and thermal drifts in this interval. We have 
calculated the rms of the residuals, which is a measure of the 
radio seeing in this interval. 

Furthermore, using these 30 min intervals of data, we 
have followed an approach using lag structure functions (see 
Hold awav et alj[T995l) to determine the phase screen speed 
and Kolmogorov exponent. In a logarithmic diagram, the 
structure function of the path length fluctuations increases 
linearly with lag, until it flattens out after a characteristic 



time. The slope of the linear increase is a measure of the Kol- 
mogorov exponent, 0/2, the time after which saturation oc- 
curs is characteristic of the speed of the phase screen, and the 
level at which saturation occurs is expected to be \fl times 
the rms path length fluctuations. We have measured these 
quantities in each 30 min interval by fitting linear functions 
to the two parts of the structure functions (Fig[3J. 

The simple rms of the data, and the fit to the saturation 
level of the lag structure function give two methods to es- 
timate The average of the ratio of these two estimates, 
after correcting for the factor of V2, is 0.995 ± 0.098, which 
is a near-perfect agreement. 

It is common to normalise the rms of the path length fluc- 
tuations to an interferometer baseline of 1 km, pointing to- 
wards the zenith. We therefore have multiplied all rms mea- 
surements with \/3/2 to convert the measurements to zenith 
values, and by (1000/230)' 3 , to scale the measurements 
from the seeing monitor's 230 m baseline to 1 km, where /3/2 
is the Kolmogorov exponent derived from the structure func- 
tion in each interval. 

The data presented here were observed in the period of 
16 September 2004 to 25 May 2006. There are only two no- 
ticeable periods where the seeing monitor was either not ob- 
serving or the measurements were below 20 /im, indicating a 
malfunction. These periods are 7 October 2004 to 19 Octo- 
ber 2004 (12.2 days); and 12 January 2005 to 14 April 2005 
(92.3 days). There are 9 shorter (one to five days) periods 
and many still shorter periods (typically a few hours) with- 
out data. The shortest outages tend to happen preferentially 
in summer and are associated with storms (the seeing moni- 
tor fails when the sky becomes opaque to the 30 GHz beacon, 
and under some situations when mains power had been lost). 
We have performed some analysis to convince ourselves that 
these short gaps in the data do not significantly affect or bias 
our conclusions below. 

The effect of the longer gaps is more difficult to esti- 
mate. The missing data essentially mean that some times of 
the year were sampled only once, and so are more prone to 
rare weather situations that may bias the characteristics of 
that period. 

4 Impact on the scheduling of mil- 
limetre observations 

4.1 Path length fluctuations as a function 
of time of day and time of year 

We have sorted the data by month and time of day (local time, 
AEST). The data observed within each month were divided 
into eight groups, each of which represents 3 h of a day in 
that particular month (00:00 to 03:00, then 03:00 to 06:00, 
and so on). The number of measurements in each of these 
groups is of the order of 220 (covering the range of 114 to 
336 measurements). Measurements of less than 20 fim have 
been deemed instrumental errors and have been ignored. The 
medians of the rms phase these groups are plotted in Fig- 
ure|4| It shows that throughout the year, the highest rms val- 
ues occur in the period between 12:00 and 15:00 AEST. The 
only exception is October, where the maximum occurs in the 
period 09:00-12:00. The plot also reveals that in any one 
month at night (between 21:00 and 06:00 AEST), the rms is 
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Figure 1: Comparison of the predicted and measured satellite phase on 23 September 2004. Top panel: The 
predicted (solid line) and measured (dots) satellite phase. The x axis shows time in UT hours and the y axis the 
phase in degrees. Bottom panel: The difference between the prediction and the measurements. The x axis shows 
time in AEST hours and the y axis the residual phase in degrees. The periodic structure in the residuals is due to 
the cycling of the air conditioning in the receiver enclosure. The relatively quiet sections before 9 h AEST and after 
17h AEST emphasize the increase in path length fluctuations during the day (cf. Fig.|4jl. 
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Figure 2: Example of seven days of seeing monitor 
data taken between 14 and 21 May 2005. The labelled 
major tic marks correspond to 00:00 h AEST. The diur- 
nal variations can be clearly seen, with relatively quiet 
nights and pronounced peaks around noon and in the 
early afternoon. Note that these values represent the 
seeing at 60° elevation on a 230 m baseline. 



lower than 650 /im, which is lower than the maximum rms in 
any month. This is interesting because it suggests that from 
the perspective of atmospheric seeing, summer nights are as 
good as, or better than, winter days. Hence it may be fea- 
sible to extend the period in which millimetre projects are 
scheduled. However, from Figure |4| it is not easy to tell by 
how much summer nights are better. We quantify this in the 
following section. 

4.2 Quantitative analysis of path length 
fluctuations 

At the ATCA, the "millimetre season", i.e., the time of year 
during which observations at millimetre wavelengths (cur- 
rently frequencies of 86 GHz and higher) are scheduled, is 
early May to August. Observing time in this period is in 
high demand, and to ensure the best possible usage of observ- 
ing time it is desirable to know to what extent the millimetre 
season can be expanded into the shoulder seasons. The see- 
ing monitor measurements from the millimetre season should 
therefore be compared to measurements taken throughout the 
rest of the year. 

Furthermore, we deem typical daytime conditions to oc- 
cur between 09:00 and 12:00, and between 15:00 and 18:00. 
The rms is remarkably worse between 12:00 and 15:00 in 
all months, and not representative of typical daytime condi- 
tions. Nighttime conditions appear to occur between 21:00 
and 06:00, where the path length fluctuations are relatively 
similar. Hence, we compare measurements taken in the mil- 
limetre season months (May to August) between 09:00 and 
12:00, and between 15:00 and 18:00, to measurements taken 
in all other months between 21:00 and 06:00. 

A measure of what fraction of time in these periods is 
useful can be obtained by estimating the amount of decor- 
relation arising from atmospheric seeing. While the amount 
of decorrelation observers are prepared to accept may vary, 
10 % appears to be a reasonable amount which would allow 



one to carry out many experiments. From expression 
the corresponding rms of the path length fluctuations fol- 
lows to 255 Lim for a wavelength of 3.49 mm (a frequency 
of 86 GHz). This needs to be scaled to a 1 km baseline to 
be comparable to our measurements, which requires a typ- 
ical baseline length used in millimetre observations, d, and 
an estimate for the Kolmogorov exponent. We assume d = 
150 m, and use the median Kolmogorov exponent of all mea- 
surements taken in the months of May to August (the mil- 
limetre season), which is f3/2 = 0.57. The rms of the path 
length fluctuations then scales to 255 /im x ( ^ ™ ) 0,57 = 
752 pm. 

For each month, we have constructed cumulative his- 
tograms of the data in each 3 h period of the day, yielding the 
fraction of measurements below a given rms value. Samples 
are reproduced in Figure [5] For example, in the top panel, 
representing all measurements of June, the "9h-12h" line in- 
tersects with the horizontal "0.6" line at an rms of 560 /im, 
which means that, statistically, in June between 09:00 and 
12:00, the path length fluctuations are smaller than 560 /mi 
during 60 % of the time. 

From the cumulative histograms we have obtained the 
fraction of measurements that lie below 750 /im. These rms 
values are listed in Table Q From May to August, between 
09:00 and 12:00 and between 15:00 and 18:00, the mean of 
the fraction of measurements which are smaller than 750 /im 
is 0.82 ± 0.05. This means that in typical conditions on win- 
ter days, during 82 % of the time the path length fluctuations 
are smaller than 750 /im and the decorrelation at 86 GHz on 
baselines of 150 m or shorter is 10 % or less. In the rest of 
the year at night, i.e., between 21:00 and 06:00, the mean is 
0.87 ± 0.07. This means that the likelihood to have see- 
ing conditions suitable for millimetre observing is a little 
higher outside the millimetre season at night, than it is in 
the millimetre season during the day. The effect is particu- 
larly pronounced in September, October, and April, whereas 
the nights of the midsummer months December and January 
are slightly worse than the average millimetre season days. 

We note that if the measurements taken between 12:00 h 
and 15:00 h in the millimetre season months were taken into 
account, then the fraction of time during which the path length 
fluctuations in this period are better than 750 /im decreased 
to 0.71 ±0.16. Then, nighttime observations outside the mil- 
limetre season appeared even more reasonable. 

4.3 Other considerations affecting the schedul- 
ing 

The scheduling of observations is not only constrained by 
weather, but is a complex process which is influenced by 
many factors. 

• Array configurations required for centimetre and mil- 
limetre observations tend to be different. Millime- 
tre observations generally require compact configu- 
rations with spacings no more than a few hundred 
metres whereas centimetre observations generally re- 
quire extended arrays. Because the compact configu- 
rations are not of interest to centimetre observers, it is 
generally more efficient to schedule the compact con- 
figurations at a time when both the days and nights 
can be used for millimetre observations. 
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RMS phase=265 microns; phase screen speed=4.2 m/s; Kolmogorov exponent beta/2=0.49 
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Figure 3: Typical lag phase structure function and a fit to these data. The data are taken on 16 June 2006. 
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Month 


00:00 h- 
03:00 h 


03:00h- 
06:00 h 


06:00 h - 
09:00 h 


09:00 h- 
12:00h 


12:00 h- 
15:00h 


15:00 h- 
18:00h 


18:00h- 
21:00h 


21:00h- 
00:00 h 


January 


0.81 


0.88 


0.79 


0.37 


0.24 


0.35 


0.73 


0.72 


February 


0.83 


0.77 


0.76 


0.38 


0.15 


0.33 


0.73 


0.72 


March 


0.91 


0.89 


0.80 


0.33 


0.22 


0.44 


0.72 


0.78 


April 


0.94 


0.93 


0.93 


0.65 


0.37 


0.76 


0.90 


0.94 


May 


0.97 


1.00 


0.99 


0.82 


0.53 


0.90 


0.99 


0.98 


June 


0.94 


0.92 


0.92 


0.78 


0.55 


0.81 


0.89 


0.92 


July 


0.92 


0.94 


0.94 


0.82 


0.56 


0.84 


0.95 


0.94 


August 


1.00 


0.98 


0.97 


0.74 


0.38 


0.82 


0.97 


0.96 


September 


0.91 


0.93 


0.84 


0.48 


0.29 


0.74 


0.90 


0.96 


October 


0.92 


0.88 


0.72 


0.31 


0.40 


0.71 


0.89 


0.93 


November 


0.89 


0.91 


0.66 


0.24 


0.16 


0.58 


0.82 


0.81 


December 


0.87 


0.88 


0.61 


0.20 


0.20 


0.42 


0.77 


0.83 



Table 1 : Fraction of time where the rms of the 
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Figure 4: Rms of the path length variations as a func- 
tion of time of day. 



• The seeing is not the only limiting factor for millime- 
tre observations. Atmospheric opacity, a strong func- 
tion of the atmosphere's total content of water vapour, 
is generally higher outside the winter, and reduces the 
sensitivity of millimetre observations. 

• The elapsed time between 2 1 :00 and 06:00 is only 9 h. 
This is not sufficient for a full synthesis observation 
unless the array is in one of the hybrid configurations, 
which are not frequently scheduled outside the mil- 
limetre season. Hence these times would not be suit- 
able for imaging observations of complex structures, 
but might be for detection experiments or monitoring 
of flux densities. 



Furthermore, the amount of data is still small, which 
makes our analysis susceptible to unusual weather condi- 
tions. For example, May in 2005 was unusually dry, and 
June to August 2005 experienced a once in 10 year wet spell, 
which may make the winter appear worse than it is on aver- 
age. 



path length fluctuations is smaller than 750 /im. 

5 Statistics of the speed of the phase 
screen 

Our analysis also yields a measure for the speed of the phase 
screen, v E . In Figure|6| we have plotted the median value of 
v s as a function of the path length fluctuations, which were 
binned into sections 50 /im wide. One can see that while the 
lowest screen speeds occur at the lowest rms values, v s is 
remarkably constant between 500 /J,m and 2400 /im, above 
which the v s increases. The interpretation of this diagram is 
that the very best rms values occur when the air above the 
observatory is extremely still. The lowest rms phase occurs 
on still winter nights, when an inversion layer has formed in 
the atmosphere: this is reflected in the low screen speeds. On 
the other hand, the highest rms phase occur during storms, 
and this seems to be reflected in high screen speeds. 

It is interesting to compare the phase screen speed with 
the physical wind speed both at the ground and at differing 
heights. We have compared the phase screen speed with the 
median ground wind speed simultaneously measured at the 
ATCA weather station. However we find no clear relation- 
ship between the two. We have also compared them in a 
statistical manner. Figure gives a cumulative histograms 
of the ATCA ground windspeed and the phase screen wind 
speed inferred from the seeing monitor. On average, the 
phase screen speed is greater than the ground wind speed. 

We are not able to directly compare the phase screen 
speed with wind levels at different heights above the ATCA. 
Instead we have obtained wind readings with altitude mea- 
sured by the Bureau of Meteorology's radiosonde program 
at the township of Moree. Moree is approximately 100 km to 
the north of the ATCA, and is in a topographically similar set- 
ting. Like the ATCA, Moree lies on the western plains, with 
the prevailing weather normally coming from the south-west 
or north-west. Radiosonde measurements from 1987 to 2005 
were used in the comparison. Analysis of the radiosonde 
data shows that the wind is normally comparatively low at 
the ground, that the wind increases significantly as soon as 
the radiosonde balloons leave the ground, but then remains 
fairly uniform with altitude up to about 2000 m above the 
ground. Comparison of the ground wind speed at Moree 
and the ATCA shows that Moree is somewhat windier place. 
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Figure 6: The speed of the phase screen as a function of 
seeing monitor rms phase. While the lowest and high- 
est screen speeds are associated with very low and high 
rms values, respectively, the screen speed is otherwise 
largely independent of rms. 
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Figure 5: Top Panel: Cumulative histogram for each 
3 h bin of the day in June. Middle Panel: Cumulative 
histogram for each bin of the day in November. Com- 
paring the two uppermost panels illustrates the gener- 
ally better seeing conditions in winter. Bottom panel: 
Comparison of typical daytime conditions in June and 
typical nighttime conditions in November. Summer 
nights appear to be slightly better than winter days in 
this diagram. 



However the Moree ground wind speeds are still lower than 
the ATCA phase screen speeds. 

Figure^also shows the cumulative histogram of the wind 
speed measured by Moree radiosondes at heights between 50 
and 200 m above the ground. This shows significantly higher 
speeds than the ground wind or the phase screen speeds. 
As the phase screen speed behaviour is more similar to the 
ground wind speed, this suggest that the phase screen is close 
to the ground. 

This conclusion is consistent with an excess amount of 
water vapour seen in the Moree radiosonde measurements: 
whereas models of water vapour content usually suggest that 
water vapour concentration falls off exponentially with a scale 
height of ~2 km, the Moree radiosonde data show an extra 
component close to the ground. 



ATCA ground wind speed - 
Phase screen speed 
Moree wind speed: (1=50-200 m - 



Figure 7: A comparison of the cumulative histograms 
of the ATCA ground wind speed, the ATCA phase 
screen speed and the speed of the winds at Moree at 
an height of 50 to 200 metres above the ground. The 
steps in some of the curves reflects quantisation in the 
measurements. The jaggedness of the Moree data is 
accentuated by comparatively modest number of mea- 
surements in this height band. 
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6 Statistics of the Kolmogorov ex- 
ponent 

Figure|8|and|9|give the Kolmogorov exponent, /3/2 as a cu- 
mulative histogram and versus rms phase. The latter figure 
gives the median binned in the same fashion as the previous 
section. The figures show that in the best seeing conditions 
/3/2 is 0.4 or even a bit lower, and so is near the value of 1/3 
expected for 2D turbulence. The exponent then rises gently 
with increasing rms phase, and reaches a value of a bit over 
0.8 in the worst seeing conditions. A Kolmogorov exponent 
of 5/6 = 0.83 is expected in the case of 3D turbulence. The 
median Kolmogorov exponent in our data is 0.56. This sug- 
gests that conditions are in a transition domain between 3D 
and 2D turbulence, and implies that the phase screen depth is 
comparable to the baseline length of 230 m. 




1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Kolmogorov exponont ((5/2) 

Figure 8: A cumulative histogram of the Kolmogorov 
exponent, (3/2. 
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Figure 9: The Kolmogorov exponent as a function of 
seeing monitor rms phase. 



7 Conclusion 

We have presented a description of the ATCA's seeing mon- 
itor along with an analysis of 21 months of data. The main 
conclusions of this analysis are that: 



• the seeing conditions on summer nights are not sub- 
stantially different from winter days; 

• the turbulence seems to be concentrated quite near the 
ground; and 

• the turbulence appears to be in transition between be- 
ing described by 2D and 3D models. This suggests 
the turbulent layer is comparable to the seeing moni- 
tor baseline length. 

The latter two points build a consistent picture of the tur- 
bulence being substantially confined to a layer in the lower 
~ 50 — 100 m of the atmosphere where the interaction with 
the ground is likely to be strongest. 

Appendix A: Decorrelation of inter- 
ferometer measurements due to at- 
mospheric seeing 

The expectation value for visibility measurements of a two- 
element interferometer, (V m ), is related to the true visibili- 
ties, V, by 

(V m ) = Ve~ a V 2 (2) 

iThompson et alJ200lL eq. 13.81), where ct,/, denotes Gaussian- 
distributed random phase fluctuations. Random phase errors 
of 30° therefore lead to a degradation in observed visibility 
amplitude to 87 % of the true value. Phase errors arising from 
water vapour scale linearly with frequency, hence to convert 
seeing monitor readings into a phase rms one can write 

°> = 27T^, (3) 

where A is the wavelength observed, and L is the rms of 
the path length fluctuations. Decorrelation is also a func- 
tion of the projected baseline length, d, and the structure of 
the atmospheric turbulence. For short baselines, the tropo- 
spheric water vapour fluctuations can be described as three- 
dimensional Kolomogorov turbulence, where the amplitude 
of the turbulence scales with the size of the turbulence cells 
as (3/2 = 5/6, such that 

^t(isib) ■ (4) 

where the 230 m is the baseline length of the seeing monitor. 
Letting V = 1 and solving for L as a function of (V m ), one 
gets 



L= — ^/-21n{V) (5) 
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